Sgr A* Companion SO-2: A Probe of Very High-Mass Star 

Formation 



Andrew Gould 

Ohio State University, Department of Astronomy, Columbus, OH 43210-1173 
gouldOastronomy . ohio-state . edu 
and 

Alice C. Quillen 

University of Rochester, Department of Astronomy, Rochester, NY 
aquillenOpas . rochester . edu 

ABSTRACT 

The star SO-2, which is orbiting Sgr A* with a 15-yr period, almost certainly 
did not form in situ. We propose that it was injected into this close orbit by 
the tidal disruption of a massive-star binary, whose primary was more massive 
than SO-2 and at least 60 M & . From numerical integrations we find that 1-2% of 
incoming binaries with closest approach equal to 130AU leave the secondary in an 
orbit with eccentricity within 0.01 of that of SO-2. If additional stars are found 
orbiting Sgr A* with relatively short periods, they could be used to probe the 
formation of massive stars in the Galactocentric region, even though the massive 
stars themselves have long since perished. 

Subject headings: binaries: general - stars: early-type - stars: kinematics - stars: 
fundamental parameters (mass) — Galaxy: center 



1. Introduction 

Almost a decade of high resolution infrared observations of the Galactic center have 
revealed at least three stars that show substantial proper-motion acceleration due to the 
supermassive black hole (BH) associated with Sgr A* (Eckart & Genzel 1996; Ghez et al. 
1998; Schodel et al. 2002; Ghez et al. 2003). The most striking of these is SO-2, which has 
a peribothron of only q = 125 AU and a period of only P = 15 yr. SO-2 has a number of 
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important potential applications. For example, by combining radial- velocity and astrometric 
observations of SO-2 one could measure the mass, the distance, and the three components 
of motion of Sgr A* (Salim & Gould 1999), which is believed to lie almost exactly at the 
Galactic center. 

However, the most remarkable characteristic of the Sgr A* /SO-2 "visual binary" is that 
it exists at all (Ghez et al. 2003). SO-2 cannot have formed in place: the gas cloud from which 
it would have condensed could never have been larger than ~ (mi/M©) 1 / 3 AU or it would 
have been tidally disrupted. Here m\ is the mass of SO-2. Hence, it must have formed at a 
much larger distance from Sgr A* and was either originally on a highly eccentric orbit that 
took it within ~ 125 AU of Sgr A* or was perturbed onto such an orbit. Then, it would have 
had to have been perturbed again at peribothron to bring its eccentricity from (l-e)<l 
to (1 — e) ~ 0.13. It is this last step that is problematic. 

Spectra at 2/im taken by Ghez et al. (2003) exhibit Br7 and Hel absorption lines and 
lack a CO bandhead. These features suggest that SO-2 is a 08-BO dwarf star and so is 
massive, mi ~ 15 M , and young, < lOMyr. The youth of the star puts a strong limit on 
the length of time available to scatter the star into its current orbit since formation. The 
dynamical friction time for such a star in an orbit of semi-major axis a ~ 10 4 AU where it 
might plausibly form, would be t^i ~ 300 Myr (a/10 4 AU)~ 1//4 , even assuming that the stellar 
profile continued as a power law down to very small radii, as suggested may be the case by 
Alexander (1999). Since several dynamical friction times would be required to achieve the 
orbit of SO-2, this appears impossible 1 . 

In this paper we propose that the orbit of SO-2 was a result of the disruption of a 
massive-star binary in the tidal field of the massive black hole at the Galactic center. When 
the binary disrupts, one star is ejected with energy greater than its incoming energy, leaving 
the other star, SO-2, in a bound orbit with the black hole. In § 2, we show analytically 
that for this mechanism to work, the former binary companion of SO-2 must have been very 
massive, of order 100 M . In § 3, we confirm this estimate by numerical simulation, showing 
that the primary star must be m 2 > 60 M Q . In § 4, we discuss how the detection of additional 
short-period (P < 100 yr) "low" -eccentricity (e < 0.95) companions to Sgr A* could be used 
to probe very high-mass star formation in the inner pc of the Galactic center. In Appendix 
A, we discuss the mechanisms by which such massive binaries might be fed to Sgr A*. 



1 If the central region is populated by a cluster of stellar-mass black holes as argued by Miralda-Escude 
& Gould (2000), then the local mass density is a factor ~ 7 lower, and hence the relaxation time is a factor 
<~ 7 longer. 
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2. Analytic Estimates 

The specific binding energy of SO-2 today is 

GM bh (l-e) 

Eb = 2q (1) 

where e = 0.87 and q = 125 AU are the eccentricity and peribothron of its orbit, and 
M hh = 3 x 10 6 M Q is the mass of Sgr A*. If SO-2 started out as part of a binary on a 
very low specific-energy (i.e., roughly parabolic) orbit with peribothron ~ q, then to become 
bound with it would have to have been ejected from the binary with a velocity 



Av = ^ = y - e) ~ 210 km s" 1 , (2) 

where v = (2GM hh /q) 1 / 2 is the velocity of the binary center of mass at peribothron. If we 
now equate Av with the internal velocity of SO-2 relative to the binary's center of mass, 



Av = Grn^^ (3) 
V a bin m tot 

where mi and m 2 are the masses of SO-2 and its former companion, m tot = m x + m 2 , and 
abin is the binary's semi-major axis, we obtain 

abin = q-r ^ • 4 

(1 - e) 2 M hh m tot 

Finally, we assume that the binary can be broken apart tidally when the accelerations are 
related by, 

2GM hh a U n Gm tot 

3 = V— — , (5) 

° 3 «bin 

where r] ~ 1 is a parameter. We then find, 

-^ = ^|(1- e) 3 M bh ~ 200 ^M , (6) 
m tot 32 

and 

a bin ^ = ^(1 - e) q ~ 4^AU. (7) 

m tot 4 

Equation (6) allows us to estimate the mass of the primary star and equation (7) the binary's 
semi-major axis. Note that in the limit, m 2 ^> mi, the left-hand sides of equations (6) and 
(7) become m 2 and a^n, respectively. On the other hand, for a roughly equal mass ratio, 
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m\ ~ m 2 ~ SOO^^M©. Such a massive SO-2 would be much brighter than is observed. 
Hence, we will generally work in the regime m 2 ^> m\. 

Equation (6) implies that the former companion of SO-2 was quite massive, with mass of 
order several hundred solar masses. However, since this equation contains the parameter 77, 
which depends on the specific internal orbital parameters of the binary, numerical simulations 
are required to determine the actual range of plausible masses for the companion. 

3. Numerical Estimates 

3.1. Single Encounters 

We conduct a series of simulations, in each of which the masses m\ and m 2 are held 
fixed, but all of the other binary parameters are allowed to vary. The three binary orbital 
orientation angles and the orbital phase are randomly chosen. The binary eccentricity is 
drawn uniformly in e^ in (but capped e^n < 0.95), and the semi-major axis is drawn uniformly 
in log a for 1 AU < a < 10 AU. We find that collisions leading to secondary orbits with the 
properties of SO-2's orbit, q = 125 AU and e = 0.87, have initial peribothra that are about 5 
AU larger, i.e., g init ~ 130 AU. All simulated collisions therefore begin with the binary on a 
parabolic orbit with g init at this value. Each set of simulations contains 20,000 integrations 
of the entire encounter between the massive-star binary and the massive black hole. 

For each set of simulations (determined by choice of mi and m 2 ) we examine the cumu- 
lative distributions of three binary input parameters (ab m , ebm, and cosibin) for the subset 
of collisions that lead to an S0-2-like orbit, namely 0.86 < e < 0.88. We also examine the 
final peribothra q of these orbits to ensure that they are consistent with the observed value 
for SO-2. Figure 1 shows the results for the 137 collisions meeting this condition for the 
simulation with m x = 15 M and m 2 = 150 M . Results for other mass combinations look 
similar. 

The first point to note from Figure 1 is that less than 1% of the 20,000 simulated 
collisions lead to SO-2 type orbits. Of course, this is partly because we have defined the 
acceptable eccentricity range extremely narrowly. However, many binaries simply do not 
break up, or if they do, they leave SO-2 in a marginally bound or unbound orbit. We return 
to this issue in § 3.2 

Of the parabolic collisions that do yield S0-2-like orbits, the distribution of binary 
eccentricities is consistent with being uniform in e^ in , which is the same as the parent dis- 
tribution of all collisions. Hence, the probability of producing an S0-2-like orbit does not 
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strongly depend on e^n- By contrast, this probability depends strongly on the binary's in- 
clination i, with prograde orbits (cosibin ~ 1) being much preferred over retrograde orbits 
(cosibin ~ — 1)- This is plausible, since it is easier to disrupt a prograde than retrograde 
binary. 

As predicted by equation (7), successful collisions have a hin ~ 4AU. There are es- 
sentially no successes for a^ n < 2AU. The decline in success toward higher semi-major 
axis seems to be less severe than the decline toward lower values. For this reason, we have 
conducted additional simulations with input distributions uniformly distributed of 1 AU < 
abin < 100 AU. We find that only about 10% of successful collisions have a^m > 10 AU, with 
the overwhelming majority of these having aun < 30 AU. 

Finally, we find that the distribution of final peribothra is sharply peaked near q ~ 
124 AU. This shows that it is unnecessary to simulate a broad range of ^ iriit because only 
those binaries with q- m i t ~ 130 AU have the potential to yield S0-2-like orbits. 

Figure 2 shows the number of successful collisions (again defined by 0.86 < e < 0.88) as 
a function of the "mass function" rrif = m\/m1 ot for 15 simulations with mi/M Q = 7.5, 15, 
and 30, and rr^/M© = 50, 75, 100, 150, and 200. Note that, in conformity with equation 
(6), the number of successful collisions is approximately a function of mj rather than of 
the two component masses separately. Also in conformity with equation (6) is the fact 
that this function peaks at mj ~ 200 M Q . There are essentially no successful collisions 
for rrif < 40 m Q which, for rrii ~ 15 M Q , corresponds to a minimum mass for the former 
companion of m 2 > 60 M & . 

As mentioned above, the absolute rate of successful collisions is quite small. Part of the 
problem is our extremely narrow definition of "success": of course the a priori probability 
of exactly reproducing the a posteriori known orbit of SO-2 is vanishingly small. One might 
plausibly argue that if SO-2 had any orbit with q = 125 AU and e < 0.9 it might well have 
been measured. Such broadened acceptance would increase the number of "successes" by 
about a factor 4 near the peak. The resulting fraction of successes is still fairly low. 

3.2. Repeated Encounters 

One should keep in mind, however, that the orbital parameters of the incident binary 
can be significantly changed even if it is not disrupted. As discussed in Appendix A, we 
expect that the incoming binary would not actually be on a parabolic orbit but rather on 
a highly elliptical one with a semi-major axis a ~ 10 4 AU, at which location a dense young 
cluster would be tidally disrupted. This induces only a slight (~ 5.4%) reduction in the 
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energy requirements of the collision. Within the formalism of § 2, the effect is equivalent to 
decreasing (1 — e) by 5.4%, or increasing e from 0.87 to 0.877. This reduces the overall mass 
scale for mj = m\/m\ ot by a factor (1 — 0.054) 3 ~ 0.85 (see eq. [6]). Given the high mass 
scale of this equation and of Figure 2, this reduction is welcome, but modest. 

However, this finite semi-major axis has another important effect: the binary re-encounters 
Sgr A* on timescales of ~ 10 3 yrs. When it does so, it preserves all the binary parameters 
it acquired in the last collision, save the internal orbital phase relative to the phase of the 
orbit around Sgr A*, which is effectively randomized. Hence, thousands of such encounters 
may take place, which increase the chance that the binary will be disrupted. Partly, the 
internal binary orbit changes each time. However, even if the orbit is left unchanged, the 
mere randomizing of the internal orbital phase permits a new opportunity for a disruptive 
encounter. 

To explore this effect, we simulate repeat encounters by binaries in orbit around Sgr A*, 
continuing the simulation until the binary is disrupted or until it completes 1000 orbits. In 
order to permit efficient computation, we place the binaries in orbits with semi-major axis 
a = 3000 AU. Then, in order to be able to compare directly with the previous simulations 
we consider as a success, collisions with the same energy change, that is, with final orbits 
having eccentricity e = 0.847. Finally, to minimize the role of Poisson fluctuations, we focus 
on the relative number successes from the entire (up to 1000 orbit) simulation with the 
number coming from the first encounter. For m\ = 15 M and m 2 = 75 M , the 1000-orbit 
simulations doubles the number of successes, while for m<i = 100 M the improvement rises 
to about a factor 2.5. 

About 90% of the post-first-orbit successes occur within the first 100 orbits, and none 
past 200 orbits, indicating that our 1000-orbit simulations were quite adequate. The great 
majority of additional successes have initial semi-major axes 1AU < <2bm < 2.5 AU. Of 
those that do not, all but one are in retrograde orbits. That is, the additional successes arise 
because binaries with unfavorable parameters are gradually perturbed into a more favorable 
regime. Given the enhanced probability for binaries with in the lowest range simulated, 
we conduct an additional simulation with m 2 = 100 M Q and 0.1 AU < a bin < 1 AU. However, 
we find very few successful collisions from this range. 

In brief, allowing for multiple encounters in repeated orbits roughly doubles the number 
of successful encounters, i.e., those leaving the secondary in an orbit similar to SO-2. 
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4. Discussion 

SO-2 is the only close companion of Sgr A* with well-determined orbital parameters and 
a high quality spectroscopic identification. However, there could be many other remnants of 
disrupted massive-star binaries still waiting to be discovered. We found from our numerical 
integrations that binaries were more often disrupted leaving a star in a bound orbit with an 
eccentricity exceeding that of SO-2. We argued in § 3 that the binary from which SO-2 was 
ejected had a mass function of at least m/ = ra\/ml ot ~ 40 M . However, from equation 
(6), binaries of this rrif should typically eject stars into orbits with eccentricity e ~ 0.92. 
For orbits with the same q = 125 AU as SO-2, this would imply a period longer by a factor 
[(1 — 0.87) / (1 — 0.924)] 3 / 2 ~ 2.2, i.e., more than 30 years. While such a star may have already 
been identified in near-infrared images, it would not yet be known whether the star would 
experience close approaches to Sgr A*. For example, the stars SI, S2 and S8 have estimated 
orbital periods of around 2000 years. However, their eccentricities are not well constrained 
(Eckart et al. 2002). Similarly, stars with the same eccentricity as SO-2 but having larger 
peribothra would also have periods too long for one to have yet determined their orbits. 
Hence, there may well be a substantial population of highly eccentric young stars orbiting 
Sgr A* that were injected into their present orbits by the disruption of massive-star binaries. 
If so, these remnants could be used to probe statistically the recent formation of very massive 
stars in the Galactocentric region, even though the massive stars themselves have long since 
perished. 

We thank Jordi Miralda-Escude and Dan Watson for stimulating discussions. Work by 
AG was supported by grant AST 02-01266 from the NSF. 

APPENDIX 

A. Feeding Delicacies to the Monster 

What mechanism might feed a substantial number of massive-star binaries to orbits 
with peribothra q init < 130 AU from which Sgr A* can "grab" an S0-2-like star? From § 3, 
only ~ 5% of binaries with semi-major axes in the decade 1 AU < a hin < 10 AU are devoured 
in this manner. Moreover, only a small fraction, perhaps O(10%), of massive stars will have 
a companion that is several times lighter orbiting with a^n in this range. Hence, perhaps 
only ~ 0.5% of massive stars passing within q ~ 130 AU will inject lighter companions into 
S0-2-like orbits. Since massive stars are fairly rare and the phase space permitted for SO-2 
progenitors is quite small, it would seem necessary that these progenitors be preferentially 
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funneled into this phase space if the collisions discussed in §§ 2 and 3 are to plausibly account 
for the existence of SO-2. 

Gerhard (2001) suggested that the Hel stars in the central parsec of the Galactic center 
were born in a 10 4 — 1O 6 M cluster that was formed more than 30pc away from Sag A* (like 
the Arches and Quintuplet clusters). The cluster sank inward to Sag A* due to dynamical 
friction and was tidally disrupted by the massive black hole, thus leaving massive stars within 
the central parsec. Consider a cluster of mass M c \ and velocity dispersion <r c i that has been 
formed outside the central pc of the Galactic center. By the virial theorem, it will have an 
effective radius R c \ ~ GM c \ja\. The cluster will sink by dynamical friction until it is tidally 
ripped apart at a separation from Sgr A*, r d , 

/MbhX 1 / 3 GM}(?M 2 / 4A / M cl \ 2 /3/ a cl \- 2 , A N 

rd~Rci(^F) bh 9 = 1.7 x 10 4 AU( -. f- ) ( , , ) . Al 

\M cl J ct 2 V1O 4 M / VeOkms" 1 / v ; 

For a star from this disrupted cluster to come within g init of Sgr A*, it must have angular 
momentum J < (2GM h i l q- mit ) 1 ^ 2 , and therefore transverse velocity v± < (2GM hh q init ) 1 ^ 2 / 'r d . 
Hence, 

v ciTC v q \M cl J V1O 4 M / GOkms- 1 1 ; 

where f C i rc = (GM/rj) 1 / 2 is the circular speed at breakup and v q = 6450 km s -1 is the velocity 
at peribothron. Hence, if the cluster is on a roughly circular orbit when it is disrupted, then 
essentially none of its stars will have a close passage to Sgr A*. 

On the other hand, if the cluster is on a roughly radial orbit, then a fraction / ~ 
1 — exp[(— qVq/r d cr c \) 2 /2] will have close passages. That is, 

/ ~ -ln(l - /) ~ 0.35 (4^f ^ (— ^_) 2 , (A3) 
J V JJ V1O 4 M / VeOkms" 1 / ' V ; 

so that for parameters that are plausible for a cluster formed in the deep gravitational well of 
the Galactic center, a substantial fraction of its members could come within g init ~ 130 AU 
following a radial-orbit disruption. In the above example, a massive-star cluster with M c \ ~ 
10 4 M Q would contain of order 10 2 100 M stars. Of these 1/3 would come within q = 130 AU 
of Sgr A* and as estimated above, of order 0.5% of these might kick out companions onto 
S0-2-like orbits. Over the lifetime of SO-2, several such clusters might form and disrupt. 
Assuming that these clusters were on radial orbits, they are a plausible source of S0-2-like 
stars. 

The final question then is whether it is plausible that such clusters will find themselves 
on radial orbits. To address this question, we first show that for a power-law density profile 
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p oc r~ v in a background Kepler potential, dynamical friction will tend to circularize orbits 
for v > 3/2 and make them more eccentric for v < 3/2. Since the maximum bound velocity 
in a Kepler potential is v max (r) = (2GM^i l /r) 1 ^ 2 , the amount of available phase-space scales 
as [fmax( r )] 3 oc r~ 3 / 2 . Hence, for v = 3/2, the phase-space density f(u,r) is a constant, fo, 
throughout the potential. The Chandrasekhar formula for the drag acceleration on a particle 
of mass m and velocity v is then a = — kv, where 

k = 4 ^"W" f dtt W f(u, r) In A — ^^rn^mf 
v Jo 3 

and where m am b is the mean mass of the ambient particles. Since In A depends only very 
weakly on radius, k is essentially independent of radius. 

Over the course of a single orbital period P, an object with specific energy E and specific 
angular momentum L will suffer mean rates of energy and angular momentum loss of 

,dE. Sdtv-a k f , 9 nl ,_, , . _ s 

= P =-pfdtv 2 = 2kE (A5) 

and 

. dL . <£ dt r x a k f , T , T ... 

(— ) = * — = --<LdtL = -2kL. (A6) 



dt' P P . 

Since 1 — e 2 = —2L 2 E/(GM hh ) 2 , the evolution of the eccentricity can be written, 

^(i-«') =2 WjO + <jgA!0 = _ at + at = a (A7) 

(it L E 

For steeper profiles, the phase-space density is higher than average at pericenter, hence the 
additional drag there tends to circularize the orbit, while for shallower profiles the higher 
phase-space density at apocenter tends to make them more eccentric. 

Over most radii, the stellar mass profile around Sgr A* is a power law with v ~ 1.8. 
Hence, there is a weak tendency to circularize orbits. If there were a core in the mass 
distribution of radius ~ 3", corresponding to ~ 2500 AU, then this core region could drive 
clusters toward eccentric orbits. That is, once the peribothron fell below the core radius, 
the declining dynamical friction inside the core (relative to that due to a power law) would 
drive the orbits to greater eccentricity and so still lower peribothron. 

However, from equation (Al), it is clear that this process could not continue very long 
before the cluster was disrupted, unless the cluster had rather extreme parameters. On the 
other hand, from equation (A2), the cluster orbit must be highly eccentric if a significant 
fractions of its disrupted contents is to reach g init . Dynamical friction by itself probably 
cannot produce such orbits. If a cluster had the required orbit to deliver SO-2 progenitors to 
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Sgr A*, it must have been born on it. This in itself is not implausible: it is certainly possible 
for clusters to form from colliding clouds of gas whose angular momenta roughly cancel. 
Whether this happens frequently or infrequently is today, however, a matter of speculation. 
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m 1 = 15 m 2 = 150 




Fig. 1. — Cumulative distributions for "successful" tidal collisions between a binary with 
masses (mi, m 2 ) = (15 M Q , 150 M Q ) and Sgr A*, with "success" defined as mi being injected 
into an orbit with eccentricity 0.86 < e < 0.88 (similar to SO-2). Shown are the semi-major 
axis Obm, the eccentricity ebm and the inclination % [n of the incoming binary, as well as the 
peribothron, q, of the final m x orbit. The simulated encounters were drawn uniformly in 
logabin, ej; in , and cosibm- Hence "success" strongly favors prograde (cosibin ~ 1) orbits, and 
semi-major axes Obi n ~ 4 AU, but does not depend strongly on e^n- All simulated encounters 
had (fc n it = 130 AU, while the final peribothra are narrowly distributed around 125 AU. 
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Fig. 2. — Number of "successful" tidal collisions (out of 20,000) between a binary with 
masses (m 1? m 2 ) and Sgr A*, with "success" defined as mi being injected into an orbit with 
eccentricity 0.86 < e < 0.88 (similar to SO-2). One may infer from eq. (6) that this success 
rate should be a function only of the combination rrif = m\jm\ Q ^ with m tot = mi + m 2 . 
The figure bears out this conjecture. The red, yellow, green, blue, and black points refer 
respectively to m 2 = 50, 75, 100, 150, and 200 M Q . In each case, trials were run for mi = 7.5, 
15, and 30 M . The success rate peaks near mj ~ 200 M , also in agreement with the naive 
prediction of eq. (6). 



